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It has a pivotal role in medical science and in industry to concentrate the acoustic energy created with 
piezoelectric transducers (PTs) into a specific area. However, previous researches seldom consider the focal 
resolution, whose focal size is much larger than one wavelength. Furthermore, there is to date no such design 
method of PTs that allows a large degree of freedom to achieve designed focal patterns. Here, an active and 
configurable planar metasurface FT prototype is proposed to manipulate the acoustic focal pattern and the 
focal resolution freely. By suitably optimized ring configurations of the active metasurface PT, we 
demonstrate the manipulation of focal patterns in acoustic far fields, such as the designed focal needle and 
multi foci. Our method is also able to manipulate and improve the cross-sectional focal resolution from 
subwavelength to the extreme case: the deep sub-diffraction-limit resolution. Via the acoustic 
Rayleigh-Sommerfeld diffraction integral (RSI) cum the binary particle swarm optimization (BPSO), the 
free manipulation of focusing properties is achieved in acoustics for the first time. Our approach may offer 
more initiatives where the strict control of acoustic high-energy areas is demanding. 



Research on acoustic focusing has led to various applications such as nondestructive testing techniques that 
inspect materials for hidden flaws^ l Usually PTs are the most commonplace devices serving as the 
actuators of acoustic focusing. As the mechanism, when an electric field is applied across piezoelectric 
materials, the polarized molecules will align with the electric field, causing the materials to change dimensions''^. 
Apart from the industrial applications, acoustic focusing utilized in medical science contributes significantly to 
therapeutic techniques as well. Ultrasound waves excited by PTs are capable of transmitting energy inside a body 
for medical purposes such as diagnostic sonography^. Other examples include focused ultrasound surgery that 
generates localized heating to treat tumors^'^, and lithotripsy that breaks up kidney stones^. 

To avoid the bulky size of a curved PT, the flat annular Fresnel PT has been invented over decades, reducing the 
volume of piezoelectric portions into a flat layer^°. However as the sacrifice, Fresnel PT intrinsically cannot 
concentrate the excited acoustic energy completely, as it is always accompanied with higher- order diffraction. 
In detail, a planar PT using a Fresnel equal-spaced array will inevitably generate the parasitic multiple divergent 
beams and higher-order convergent beams, making the focal spots less applicable. Actually, there has been to date 
no such acoustic design technique that allows us to achieve arbitrarily designed focal pattern along the axis. For 
example, one expected focal pattern for ultrasonic surgery is a specific segment of high acoustic energy along the 
axis (both its distance away from the PT and its focal depth are designable), i.e., an acoustic far-field focal needle, 
which was never obtained by PTs. A finite-length focal-needle pattern is also quite promising for particle 
operation and acceleration, which was developed in optics^^'^^. 

Additionally, the focal resolution created with traditional PTs is usually low, whose focal size is much larger 
than one wavelength L Since the focal resolution can be improved with the wave frequency increased, previous 
researches usually ignore the consideration of improving the resultant focal resolution. However, it is noteworthy 
that an excitation of a higher frequency demands more energy consumption and suffers stronger attenuation. 
Besides, acoustic aberration could also severely blur the focal resolution. Thus, the rational improvement is to 
increase the relative focal resolution with the same excited frequency remaining. 

In this paper, we establish the acoustic PT prototype, whose piezoelectric elements are squeezed into a flat thin 
layer compared to the scale of the entire device. The active planar interface also extends the knowledge of acoustic 
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metasurface engineering for the deflection of sound beams using 
passive elements^^'^^. Through the optimized ring configurations of 
the active metasurface PT, we are able to manipulate the focal pattern 
and the focal resolution in acoustic far fields. Firstly, we design the 
far-field finite-length focal needle with the manipulated distance and 
depth. Its focal resolution is subwavelength for the full width at half 
maximum (FWHM), and it propagates without divergence for a 
distance of —5.9/1 as designed, longer than the depth ~ 4/1 of the 
reported optical needle^^. These two designed focusing properties 
created with PTs were never achieved in acoustics, to the best of 
our knowledge. To further verify the robustness of our manipulation 
of the focal pattern, via another optimized ring configuration, we 
obtain the designed far-field multi foci, whose FWHM (~0.45>^) 
beats the Rayleigh diffraction limit of conventional acoustic instru- 
ments (0.5/1). Besides, to demonstrate the manipulation of the focal 
resolution, we design the extreme case of the super- oscillatory super 
resolution, whose size is —0.3/1 in acoustic far fields, much smaller 
than the diffraction limit. 

In the three-dimensional view of the configurable planar metasur- 
face PT prototype in Fig. 1(a), piezoelectric rings (red) are unevenly 
spaced with hard-boundary mask rings (blue) in between. A type of 
common artificial ceramic is employed as the piezoelectric material: 
lead zirconate titanate PZT-5H^^. In the radial view, the thickness q is 
set identical for all PZT-5H rings, and the ring configuration (rl, rrl, 
r2, rr2 ...) will be optimized according to different focusing manip- 
ulation. The thin hard-boundary mask rings in between, through 
which no sound can pass, are the complements of the spaced gaps 
between PZT-5H rings, co-planar with z = 0. The entire PT is axi- 
symmetric with respect to +z, and the upper surface of the structure 
at z = 0 can be regarded as a flat active metasurface according to the 
radial cross -sectional view. In our following demonstrations in air 
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Figure 1 | (a) Schematics of the configurable planar metasurface PT 
prototype. Inside the dashed box is the radial cross-sectional view of the 
ring configuration, showing the unevenly- distributed piezoelectric 
elements and the hard boundaries, (b) The B.C.s of each piezoelectric ring 
observed from the radial cross-sectional view. 



(density: po = 1.21 kg/m^; speed of sound: Cq = 343 m/s), we will 
show the designed focal pattern and the focal resolution created with 
the PT prototype in acoustic far fields, the simulation of which is 
carried out by the finite element method (FEM). In detail, COMSOL 
Multiphysics is the platform we use, and the simulation is facilitated 
by the coupling of the embedded acoustic module and the acoustic- 
piezoelectric interaction module concurrently. 

The design method and the related physics are elaborated in this 
part. The electromechanical constitutive equations governing the 
piezoelectric effect of PZT-5H are written in the stress-charge form^^: 



T = CeS- 



D = eS-\- £sE, 



(1) 



where T and S are the vectors of stress and strain; E = — is the 
electric field, which is rephrased with the electric potential (j) exerted 
on the PZT-5H rings; Ce the elasticity matrix evaluated at a constant 
electric field; e the coupling matrix; D the electric displacement; Ss the 
electric permittivity matrix evaluated at a constant mechanical strain. 
Usually, Ce is straightforwardly given as an anisotropic symmetric 
matrix; 8s = SoSrs includes the relative permittivity matrix £^5. These 
parameters of PZT-5H are listed in Supplementary Information. For 
each PZT-5H ring whose poling direction is aligned toward +z, the 
boundary conditions (B.C.s) are indicated in Fig. 1(b). The structural 
B.C.s for the inner and the outer sides are free of constraint, while the 
bottom is fixed with the structural displacement u = 0. The top 
undergoes the interaction between sounds and structures'^: 



T-ns 



Po 



(2) 



S=-[(Vu)^ + Vu], 



where T and S are the tensors of stress and strain; ris is the outward- 
pointing unit normal vector seen from inside PZT-5H. Briefly, the 
first equality indicates the acoustic pressure loadp; the second makes 
the normal acceleration of p on the boundary z = 0 equal to that of 
the structural displacement; the third is the intrinsic constitutive 

equation of S and u inside solids. Meanwhile, the electrical B.C.s 
for the inner and the outer sides are free of charge • D = 0, whilst 
the bottom touches ground (j) = 0. The top is assigned with mono- 
chromatic (j) = Vocos(2nfi). One common technique to impose the 
voltage bias is a sandwiched structure with PTH-5H elements 
between two Au electrode layers. 

Next, in order to manipulate the far-field focal pattern or the focal 
resolution created with the planar metasurface PT prototype, we 
propose the acoustic RSI in conjunction with BPSO'^"^°. The deriva- 
tion starts from the Kirchhoff-Helmholtz integraP^: 



[p(r',a;)V'^(|r-r'|)-^(|r-r'|)V'p(r',co)]-n'dS'=p(r,a;), (3) 



where r e Q, ^(x) = exp{ikox)/{4nx) and/co = 2nflcQ.p{') and^(-) are 
monochromatic wave functions defined throughout the domain Q 
bounded by a closed surface S, while n' • V denotes the differentiation 
along the inward normal to S. Provided that Q is the half space z > 0, 
S thus consists of the plane z = 0 and a hemisphere of infinite radius, 
centered at the origin. After adopting the Sommerfeld radiation con- 
dition lim [dp{Y,(jj) / dr — ikQp{Y,(jj)\ =0 to Eq. (3), we obtain: 



a exp(z/coi^i) exp(//coi^i)5p(r',ft>) . . 

.<^)^ —, — ]dx dy =p{r,a)), (4j 



'dz' 



dz' 



where Ri = ^ {x-x'f + {y-y'f + {z-z'f ,z>^. On the other 
hand, if r situates in the lower half space z < 0, we will obtain^^: 
d expiikoRi) expiikoRi) dpir' ,w) , 

[P(>* '^)^ n n ^Zr-]^^^y =0,(5) 
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Figure 2 | (a) The normalized squared absolute pressure, displaying the pattern of the designed finite-length far-field focal needle, (b) The field 
distribution of the squared absolute pressure around the focal needle. (c,d) The radial distributions of the squared absolute pressure at the cross sections z 
= 201 and z = 24/i, with their respective field distributions. 



where R2 = \J {x — x'Y + (y — y')^ + {z + z'Y ,z < 0. Subsequently, 
making use of the connection between Ri and at 2: = 0, we incorp- 
orate Eq. (4) and Eq. (5), bringing about the acoustic RSI: 



27rJ, 



Applying the first-order Born approximation^^ to Eq. (6), which 
simplifies the vibration at the fluid-solid interfaces at 2: = 0 to be 
uniform and binary, we finalize the equation as the acoustic RSI cum 
BPSO: 



\pirM\' = 



puni(r\cD) 



z exp(ikoRi) 



iko- 



dx'dy' 



,{7) 



where punii') is a binary function describing the ring configuration at 
z = 0, optimized for the designed focal pattern or focal resolution 
\p(r,w)\^. Other than the requirement of acoustic pressure as well as 
its normal derivative in Eq. (3), Eq. (7) only requires the acoustic 
pressure on the surface, suitably for BPSO. 

In general, BPSO is a nature-inspired evolutionary algorithm for 
stochastic optimization^^. At first, the designed energy distribution 
|p(r,CL)) 1^ will be preset to benchmark acoustic focal patterns or reso- 
lutions in respective scenarios. Simultaneously, by means of Eq. (7), 
BPSO will be implemented to optimize the parameter r = 
[ri,r2,r3,. . .,rAr] in order to fulfill the benchmark \p{r,w)\^. For simpli- 
city, the ring width Ar of the planar metasurface PT is fixed in our 
design. Once r is given, one can obtain another parameter rr = 
[rri,rr2,rr3,...,rr^] by rr = r + Ar. Then, the parameters r and rr 
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Figure 3 | (a) The normalized squared absolute pressure, displaying the pattern of the designed far-field multi foci, (b) The field distribution of the 
squared absolute pressure around the multi foci. (c,d) The radial distributions of the squared absolute pressure at the cross sections z = 18.4/1 and 
z = 24.2/1, and their respective field distributions. 



determine the ring configuration in Fig. 1(a). Although the fixed Ar 
sets a constraint in designing the active metasurface PT, this low- cost 
approach always works well in many applications^ ^ as we will show 
later. 

To vindicate the proposed method, i.e. the acoustic RSI cum 
BPSO, in the manipulation of acoustic focusing, we first demonstrate 
the manipulation of the acoustic focal patterns such as the designed 
focal needle and the designed multi foci. The arbitrary design of a 
focal pattern is impossible if we simply resort to acoustic wavefront 
construction by the method of effective medium^^'^^. In our case for 
the pattern of a finite-length focal needle on axis, we conveniently 
select Vq ^ 5yand/= 100 /cHz that generates acoustic waves of >^ = 
3.43 mm in the space 2: > 0. Note that for the purpose of a finite- 
length focal-needle pattern in the far field, we require a depth of 



continuous acoustic focal energy along the axis with the low energy 
level at the rest, whilst the location, i.e., both the depth of the needle 
and the specific positions of the two ends away from the transducer, 
could be subtly designed as well. In Fig. 2(a), the on-axis focal-needle 
pattern, whose position is preset to extend from 19.2>^ to 25. 1>^, is 
designed as the orange dashed curve |p(r,co)|^, while simultaneously 
the optimized ring configuration is calculated by the acoustic RSI 
cum BPSO as described above. The optimized ring configuration 
listed in Supplementary Information includes 30 PZT-5H rings with 
the maximum radius —180 mm, while q is adjusted optimally to be 
4 mm. 

In the same plot, the full -wave acoustic field generated by the 
planar and active metasurface PT as the real case is simulated by 
the FEM as the blue curve, using Eqs. (1,2), and the result coincides 
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Figure 4 | The radial distribution of the normalized squared absolute 
pressure at the cross section z = 20.06>2 on the left, showing the focal size 
of the far-field acoustic super-oscillatory super resolution. On the right is 
the corresponding field distribution of the squared absolute pressure. 

with the designed finite-length focal-needle pattern (orange dashed 
curve). It is noteworthy that the standing- wave -like oscillation from 
OX to ~10>^ along z agrees with the classic characteristic of acoustic 
near fields^^. Besides, we notice that most energy is focused into the 
designed far-field focal region from 19.21/1 to 25.1/1, implying the 
focal needle of the depth 5.89/1 as designed, longer than the depth 
~ 4/1 of the reported optical needle^^ In Fig. 2(b), the field pattern of 
|p(r,co) 1^ simulated by FEM around the focal needle is displayed. The 
contrast between the intense needle and the ambient quiet field is 
apparent, which meets the requirement of the focal-needle shape. 

It is necessary to mention the subwavelength focal resolution of 
the focal-needle pattern in Fig. 2(c,d). The radial distributions of the 
squared absolute pressure are plotted on the left while their field 
distributions are on the right. The FWHM of the focal size in the 
far field is measured to be ~0.64>^ dXz = 20/1 and z = 241, smaller 
than one wavelength but larger than the diffraction limit. We may 
also conclude from Fig. 2(c,d) that the acoustic needle pattern is 
formed of a nearly constant subwavelength width throughout. 
Moreover, the intensity of the sidelobes in Fig. 2(c,d) is drastically 
smaller than the central intensity, crucially for potential imaging 
applications. Note that such subwavelength acoustic focusing is gen- 
erated in the true far field without resorting to evanescent acoustic 
waves limited to the near field only. We believe this is the first time 
that the pattern of an acoustic finite-length focal needle is actualized. 

To further show the robustness of the acoustic-focusing manip- 
ulation, we take the example of the acoustic far-field multi foci as 
another designed focal pattern. Here, the multi foci are designed as 
the four discrete foci along the axis in the far field. The corresponding 
normalized energy pattern |p(r,(X>)|^ designed by the acoustic RSI 
cum BPSO is the orange dashed curve in Fig. 3(a), using the ring 
configuration which is simultaneously optimized in this case. Note 
that the ring configuration here is designed and optimized in the 
same way except for a different focal pattern (benchmark). It includes 
28 PZT-5H rings, whose parameters are hsted in Supplementary 
Information, q = 3 mm is optimized here while Vq and /remain 
the same. The blue curve indicates the full-wave simulation by FEM 
using Eqs. (1,2). Again, the satisfactory agreement between these two 
outcomes confirms our pattern design. The corresponding field dis- 
tribution in Fig. 3(b) is simulated with respect to |p(r,(L>) |^ around the 
multi foci. Also, we notice that the focal resolution (FWHM —0.45/1) 
of the multi foci in Fig. 3(c,d) is subwavelength and even beats the 
Rayleigh diffraction limit of 0.5>1, which was never realized in terms 
of FT technology. 

As mentioned in the above two cases, by means of the properly- 
optimized ring configurations, we are able to achieve the designed 
acoustic focal patterns with a subwavelength resolution as well as a 
sub -diffraction -limit resolution. The underlying physics of achieving 



the focal resolution beyond the restriction of diffraction limit by 
means of the multi-ring active metasurface transducer is the inter- 
ference and the diffraction of excited sound beams, which has been 
thoroughly elaborated as a wave phenomenon in the scenario of 
optics^ \ In essence, the superposition of the excited zero- order 
Bessel patterns at the focal plane breaches the traditional Rayleigh 
diffraction limit, while the algorithm BPSO simultaneously opti- 
mizes the ring configuration to further improve the acoustic focal 
resolution^ \ If we intend to further enhance the focal resolution to 
the deep sub -diffraction -limit focal size substantially smaller than 
0.38/1, the situation will turn into the extreme case of the super- 
oscillatory super resolution^^'^^. 

Our design method, the acoustic RSI cum BPSO in terms of PT 
technology, is able to manipulate the acoustic focal resolution as well 
and to increase it to the acoustic super- oscillatory super resolution. 
In this case, we choose Vq ^ 100 y and / = 40 kHz for the electric 
potential exerted on all PZT-5H rings, which produce acoustic 
waves of = 8.575 mm in air. After setting the acoustic super-oscil- 
latory super resolution (FWHM = 03k dX z = 20.06X) as the 
designed focal resolution, we simultaneously optimize the ring con- 
figuration of the active metasurface PT in the same say, which 
includes 19 PZT-5H rings as shown in Supplementary Informa- 
tion, while q = 2 mm is adopted here optimally. The optimized result 
is the orange dashed curve in Fig. 4, showing the designed radial 
distribution of the normalized |p(r,CL)) |^ dXz = 20.06L The blue curve 
shows the corresponding full-wave simulation. The field distribution 
of the squared absolute pressure at the same cross section is plotted 
on the right in Fig. 4 to exhibit the result of the focal- resolution 
manipulation. 

We manage to control the focal size of the super-oscillatory super 
resolution to be 0.3/1 in acoustic far fields, substantially smaller than 
0.5/1. Nevertheless, if the FWHM of the focal spot is smaller than 
0.38>1, there is an inevitable tradeoff uniquely in terms of the super- 
oscillatory case, which does not exist for a larger FWHM. Namely, 
the cost of the super-oscillatory super resolution is that most on-axis 
focal energy would be unavoidably squeezed away into the sidelobes 
as shown in Fig. 4^^'^\ In principle, the deep sub -diffraction -limit 
super resolution can be manipulated to be infinitesimal, but the 
sacrifice is the increasingly- outsized sidelobes that make the super- 
oscillatory super resolution less efficient. However in spite of the 
naturally-inevitable tradeoff in the super-oscillatory case, people also 
try applying it in biological imaging^^. 

Physically, the proposed active and configurable planar metasur- 
face PT used for the manipulation of acoustic focusing is a sort of 
binary- amplitude exciter for directly focusing the excited sounds into 
a certain pattern or a designed resolution in the post-evanescent 
fields, by optimally modifying the constructive interference of a large 
amount of acoustic beams. The acoustic RSI cum BPSO is applied to 
optimize the ring configuration on the active metasurface after cer- 
tain focusing targets such as focal patterns or focal resolutions are 
chosen. Additionally, the manipulation of acoustic focusing is not the 
translation from optics, whose binary mask is passive without energy 
feeding^^. Here, we fundamentally change the passive binary-ampli- 
tude baffle in optics into the active binary- amplitude exciter in acous- 
tics that transforms itself into the feeding. Furthermore, we 
demonstrate the versatility of our design method, such as the real- 
ization of the focal-needle pattern, the multi-foci pattern and the 
manipulated super- oscillatory super resolution in acoustic far fields. 
Consequently, this approach to manipulate focusing properties is 
unique in acoustics. 

In summary, we propose the method to manipulate the acoustic 
far-field focusing created with the active and configurable planar 
metasurface PT. The designed focal patterns and focal resolutions 
are demonstrated respectively. By means of the acoustic RSI cum 
BPSO, a large degree of freedom is obtained to manipulate various 
far-field focusing phenomena via optimized ring configurations. 
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Here for the first time, the pattern of a finite-length focal needle of a 
subwavelength resolution is designed and realized; another example 
is the creation of the multi-foci pattern of a sub-diffraction-limit 
resolution; for the focal-resolution manipulation, the extremely-high 
resolution, i.e., the acoustic super-oscillatory super resolution, is 
demonstrated as well using the active and configurable planar meta- 
surface PT prototype. Due to these significant advantages and ver- 
satility, the manipulation of acoustic focusing may offer a 
reformative framework in medical and industrial technology where 
the strict control of acoustic high-energy areas is demanding (bio- 
medical actuator, focused ultrasound surgery, lithotripsy, nondes- 
tructive testing, etc.). 
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